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Oh' 

\^ I We determine the time evolution of the colour screening mass in high energy 

■ nuclear collisions, as provided by the parton cascade model. Using our result, 

we discuss the onset of deconfinement and the onset of quarkonium suppression 
^ ] in a general, not necessarily equilibrated environment of strongly interacting con- 

stituents. 
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Abstract 



in 
a^ 
o 

. Statistical QCD predicts that with increasing temperature, strongly interacting media 

will undergo a transition from hadronic matter to quark-gluon plasma. The aim of high 
CLi! energy nuclear collisions is to study this transition and the resulting new deconfined 

Q_i' phase in the laboratory. However, statistical QCD deals with large systems in thermal 

^ ■ equilibrium, while nuclear collisions provide small, rapidly evolving systems. A study of 

the onset of deconfinement in a microscopic, non-equilibrium space-time picture would 
therefore be of great help in understanding if and to what extent theory and experiment 
^ I can be expected to meet. The aim of the parton cascade model [Q, |^ is to describe 

' the evolution of a high energy nucleus-nucleus collision through perturbative partonic 

interactions embedded in a relativistic transport theory. Hence it seems to be a good 
starting point for a 'dynamic study of deconfinement'. 

The initial state of two colliding nucleons can be pictured as colliding beams of confined 
partons. The incident partons have a distribution in intrinsic transverse momentum, 
leading to some average value which in turn defines an average transverse parton 

size tt = l/ikT)- When two nuclei collide at high energy, the transverse parton density 
increases considerably, since now we have a superposition of nucleon-nucleon collisions. 
In this medium, colour screening will destroy the association of partons to particular 
hadrons, since for a sufficiently high density of colour charges, the colour screening radius 
becomes much smaller than the typical hadronic scale. Hence we expect the onset of 
deconfinement for some characteristic density or for an equivalent screening scale. In 
perturbative QCD, the method to calculate the colour screening mass /i (the inverse 
of the screening radius td) is well-known; an extension to the non-equilibrium medium 
provided by cascade models has also been proposed 0, ^. We shall therefore use this 
method to obtain the time evolution of the screening mass in the parton cascade 
model. 
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Given we have to specify for what value deconfinement sets in; the parton cascade 
model itself does not identify such a point. It just has the primary collisions between the 
'confined' partons of the incident nucleons, followed by successive interactions between 
primary as well as produced partons and eventually by the hadronisation of partons ac- 
cording to some assumed scheme. There are two possible ways to specify deconfinement 
in such a picture. One can obtain at any given time the transverse density profile of the 
parton distribution, and use the percolation point of partonic discs in the transverse plane 
to define the onset of deconfinement Here we shall follow a somewhat different 

path, taking as critical screening mass fic the value obtained from lattice studies. These 
provide the temperature dependence of fi and in particular also its value at the deconfine- 
ment point Tc, where one has fic = fJ'{Tc) — 0.4 — 0.6 GeV p|. It is not obvious that such 
an equilibrium value is really applicable to the non-equilibrium situation provided by the 
parton cascade model. A good check would be to compare in this model the transverse 
parton density at fic to the percolation value. Work on this is in progress. 

Let us briefly recall the most important features of the parton cascade model. 

• The initial nucleus-nucleus system is treated as two colliding clouds of partons, 
whose distributions are fixed by the nucleonic parton distribution functions de- 
termined in deep inelastic lepton-nucleon scattering, and by the nucleon density 
distributions in the nuclei. 

• The parton cascade development starts when the initial parton clouds interpene- 
trate, and it follows their space-time development due to interactions. The model 
includes multiple elastic and inelastic interactions described as sequences of 2 — >■ 2 
scatterings, 1 — * 2 emissions, and 2 — > 1 fusions. It moreover explicitly accounts 
for the individual time scale of each parton-parton collision, the formation time of 
the parton radiation, the effective suppression of radiation from virtual partons due 
to an enhanced absorption probability by others in regions of dense phase space 
occupation, and the effect of soft gluon interference in low energy gluon emission. 

• Finally, the hadronization in terms of a parton coalescence to colour neutral clusters 
is described as a local statistical process that depends on the spatial separation and 
colour of the nearest-neighbour partons 0. These pre-hadronic clusters then decay 
to form hadrons. 

Since we are interested in the early state of the collision when the strongly interacting 
matter is still in partonic form, the hadronisation part of the model will not play a role 
for our considerations. 

To specify the further basis for our considerations, we add some details of the parton- 
parton scattering considered in the parton cascade model The elementary parton 
scatterings a + b — > c + d and the fusion processes a + b —>■ c* themselves are devided 
into two distinct classes: 

i) hard parton collisions with a sufficiently large momentum transfer or invariant 
mass s to apply perturbative QCD; and 

ii) soft parton collisions with low momentum transfer or invariant mass s, which 
will be modelled phenomenologically. 



2 



This division is necessary to regulate the colhsion integrals which appear in the transport 
equations describing the evolution of the partonic system and which are singular for small 
momentum transfer in the Born approximation. The parton-parton cross-section is 
rendered finite by invoking an invariant hard-soft division scale such that the collisions 
occuring at a momentum scale > are treated perturbatively, whereas a soft, non- 
perturbative interaction is assumed for those with < Pq- The total cross-section for 
collisions between two partons is then written as 



CTablS 



where the symbols have their usual meaning. The specific value of po is fixed by demand- 
ing that the energy dependence of the total cross-section for pp collisions is correctly 
reproduced when the expression is convoluted with the structure function of the nucleons 




within an eikonal approximation [Q]. One way to satisfy this requirement is a scale |10 



Po=Po(v^)=0.5(v/^r^ (2) 

and we shall use this form, where y/s is the nucleon-nucleon c.m.s. energy. 

In the two-term Ansatz (1), the hard and soft contributions now have to be specified. 
For the hard collisions above po we use the standard form 

-dQ^ = Y6^2\M ab-.Us,Q )| - (3) 

where |Mp is the process-dependent spin- and colour-averaged squared matrix element 
in Born approximation. To account for higher order contributions, we modify the cor- 



responding cross sections by a i^-factor of 2.5 |TT|]. Soft collisions between two partons 
are assumed to proceed through a very low-energy double gluon exchange. This provides 
a natural continuation to the harder collisions above po? where the dominant one-gluon 
exchange processes in gg —>■ gg, gq — > gq, and qq —>■ qq have the same overall struc- 
ture |T2[. A simple and physically plausible form for the soft cross-section continues the 



hard cross-section for below p^ down to = by introducing a regulating term 
so that we have 

dKb^-lcd 27ra^(pg) 

rfQ2 ~ Q2 + ^2 y ) 

Thus (3 acts as a phenomelogical parameter which governs the overall magnitude of the 
integral cr*"-^* ~ In [{pi + (3'^)/ (3\ it is estimated to be in the range of 0.3 - 1.0 GeV, and 
we shall use (3 = 0.5 GeV 0. 

The only further parameter to be specified is the virtuality cut-off for the evolution 
of the time-like partons, below which they are assumed not to radiate; this is taken to be 
1.0 GeV, as determined from fits to the hadron production in e^e~ collisions 0]. Finally 
we add that the maximum possible longitudinal spread for the gluons, Az ~ 1/xP, with 
P denoting the longitudinal momentum of the nucleon, is taken as 1 fm |ir 
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The parton cascade model provides the phase space distribution of the partons. With 
this given, we have the general form for the colour screening mass in the one loop approx- 
imation 1^, [14 1 



lim 

k1-o 



q- k 



(5) 



where is the strong coupling constant, the fi specify the phase space density of gluons, 
quarks, and anti-quarks and q runs over the flavour of quarks. It is easy to verify that in 
the case of ideal gas of massless partons, where the fi reduce to Bose-Einstein or Fermi- 
Dirac distributions (with vanishing baryochemical potential fis), Eq. (5) reduces at high 
temperatures to 



AnasT^ ( 1 + 



6 



(6) 



On the other hand, for the same system at large and low temperature, we get 



vr 



(7) 



Following 0, we shall assume that Eq. (5) holds also for the non-equilibrated partons 
created by partonic scattering and radiation in the early stages of the nuclear collision. 
We use the parton cascade model to estimate the phase-space density of the partons. The 
phase space density of the partons thus produced can be written as 



2(27r) 



gzV \k 



■fi{kT,y), 



where gi gives the relevant degeneracy with gg = 16, and gg = g-g = 6, V is the volume 
occupied by the partons at the time r, and k^ and y are the transverse momenta and 
the rapidity of the parton under consideration. The partonic distribution will be initially 
anisotropic with respect to the beam axis and thus the screening mass of a gluon in such 
a matter will depend on its direction of propagation. 

It is known that the initial distribution of the partons in the parton cascade model 
resembles a plateau extending to rapidities ±Y. Y is found to be of the order of 2 at 
^/s = 20 A-GeV, about 3 at y/s = 200 A-GeV, and more than 4-5 at energies likely 
to be attained at LHC To allow a comparison of the evolving scenarios as the energy 
is increased, we shall keep y = 2 in the following. The step-wise development of the 
parton-cascade in time as implemented in Monte-Carlo simulation VNI is of great 
help in accounting for partons at some given instant of time, in a given volume. It may be 
noted that in the parton cascade model, the partons are given a formation-time and are 
counted as real particles only after this time. For ease of computation, we approximate 
the parton distributions as 



fiikT, y) = ^U{kT) [e{y + Y)- e{y - Y)] 



(9) 



where fiikr) = dNi/dk^ for the partons in the volume we choose for estimating the 
screening mass. We are interested in seeing the variation with time r and with distance 
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tt in a plane transverse to the collision axis; in central collisions, the density of partons 
is highest along this axis. We shall therefore choose different zones Ri < vt < Ri+i in the 
central slice near z = 0. The screening masses in central A A collisions are then given by 

Bi 



A8as ( sin ^(tanhF)\ 



-/p(M + -E{/.(M + MM} 

9g Qq f 



(10) 

^" = (^ + sinhFsin-i(tanhr))- ^^^^ 

for the transverse and longitudinal directions, respectively. For Y = 2, the factor in 
brackets in Eq. (11) is about 1.2, so that fix — within 10 %, and hence we shall in the 
remainder consider only fix = f^- 

The evolution time slice of interest is somewhat arbitrary; it can probably be fixed 
more precisely when discussing specific signatures, and we shall come back to this point 
a little further on. The calculations to be shown here will start at proper time r = 0.1 
fm, assuming that at least this much time is needed to establish any kind of medium. 
For r > 0.5 fm, the screening masses tend to become only weakly time-dependent, so we 
present calculations up to about 1 fm. 

We begin by calculating the screening mass for partonic matter in central S — S 
collisions at SPS, RHIC, and LHC energies, for transverse distances < 2 fm (Fig. 1). 
Here and in all subsequent calculations, the cascade starts when the two incoming nuclei 
are centered ai z = ±1 fm and r = — 1 fm. At SPS energy, the screening mass is found 
to be always less than 0.25 GeV and thus considerably below the deconfinement point 
/ic — 0.5 GeV. We therefore expect that the partonic matter produced in this case could 
not have been deconfined. In contrast, RHIC energies bring /i just into the deconfinement 
zone, and for LHC, n » fic for all times considered here, so that at the LHC, all collisions 
at reasonable impact parameters will lead deeply into the deconfinement regime. 

Next we turn to Pb — Pb collisions at the SPS, see Fig. 2. Here there is a hot center, 
i.e., at early times, a region for sufficiently small reaches screening masses around fic- 
For Au — Au collisions at RHIC (Fig. 3), /i increases by more than a factor three, so that 
now essentially the entire collision volume falls into the deconfined region. To study the 
onset of deconfinement at RHIC would thus seem to require smaller nuclei, and so we 
show in Fig. 4 corresponding results for Cu — Cu collisions, confirming this expectation. 

Before considering applications of these results to the study of deconfinement signals, 
we should note that our calculations are based on perturbative partonic interactions and 
a phenomenological extension thereof into the soft regime. There are indications 
that non-perturbative corrections may lead to considerably larger screening masses, and 
it would clearly be of great interest to obtain more precise results from finite temperature 
lattice QCD. 

For a test of deconfinement in nuclear collisions, the production and suppression of 
quarkonium states appears so far to be the most suitable probe []T^. The tight binding 
of most of these states prevents their dissociation in confined matter [|1^ ; in a deconfined 
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medium, a hierarchy of suppression is predicted, governed by the size or the binding energy 
of the state in question ||19|. In Table 1, we summarize the screening masses obtained in 
potential theory based on a screened confining potential p!9| . 



State 




Xc 


T 




Xb 


/xf ^ (GeV) 


0.70 


0.34 


1.57 


0.67 


0.56 



Table 1: Critical screening masses for quarkonia. 

From Fig. 2 we see that at the SPS, collisions of sufficiently heavy nuclei reach a regime 
in which much of the Xc production should be suppressed by deconfinement. The effect on 
direct J/ ip production seems to be a more quantitative question - slight non-perturbative 
contributions could perhaps enhance enough to reach deconfinement values. In Au — Au 
collisions at RHIC, essentially all charmonium production should be suppressed, see Fig. 3, 
and if sufficient statistics become possible, the onset of bottonium suppression could be 
studied here for the first time. The LHC, finally, would seem to be an ideal tool to carry 
out a systematic study of bottonium suppression. 

We thus find that the effect of screening mass considerations in the parton cascade 
model leads to a pattern of quarkonium suppression which is very similar to what has been 



obtained in other, more global approaches [20|. In particular, recent suppression studies 



based on string p or parton percolation P, agree very nicely with the present results. 
An open question which could be addressed in the present, time-dependent approach is 
the dissociation of a charmonium state in an evolving medium. As a first guess, we have 
assumed that such states are dissociated once the screening length reaches deconfinement 
values. For more quantitative estimates, one should solve the bound-state problem in the 
presence of time-dependent screening, to see how long the nascent charmonium state has 
to spend in a deconfining medium before it is fully dissociated. 

In closing we want to emphasize that our main aim was to outline how a microscopic 
evolution scheme based on parton interactions can be used to study the onset of deconfine- 
ment and its effect on deconfinement signatures. We have here used the parton cascade 
model; it would clearly be of interest to see if other, similar or not so similar models 
corroborate our results. 
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Figure 1: The time evolution of the colour screening mass in central S — S collisions at 
SPS, RHIC and LHC energies. 
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Figure 2: The time evolution of the colour screening mass in central Pb — Pb collisions at 
SPS energy in different radial zones. 



9 



Screening Mass From PCM; Au+Au@ 200 A-GeV 

2.0 I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



1.5 - 



> 
CD 



1.0 - 



0.5 - 



0.0 




1 — r 



|Y|<2 



0<rT<2 fm 

2<rT<4 fm 
4<rT<6 fm 

6<rT<8 fm 



J I I I I I I I I I I I I I I I L 



0.0 0.2 0.4 0.6 0.8 1.0 

T (fm/c) 



Figure 3: The time evolution of the colour screening mass in central Au — Au collisions 
at RHIC energy in different radial zones. 



10 



Screening Mass From PCM; Cu+Cu@ 200 A-GeV 
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Figure 4: The time evolution of the colour screening mass in central Cu — Cu collisions 
at RHIC energy in different radial zones. 
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